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ABSTRACT 


Multistatic  radars  utilize  multiple  transmitter  and  receiver  sites  to  provide  several 
different  monostatic  and  bistatic  channels  of  observation.  Multistatic  passive  and  active  radar 
systems  can  offer  many  advantages  in  terms  of  coverage  and  accuracy  in  the  estimation  of 
target  signal  parameters.  Unfortunately,  their  performances  are  heavily  sensitive  to  the 
position  of  receivers  and  transmitters  with  respect  to  the  position  of  the  target.  It  is  well  know 
that  geometry  factors  play  an  important  role  in  the  shape  of  the  ambiguity  function  (AF) 
which  is  often  used  to  measure  the  possible  global  resolution  and  large  error  properties  of 
the  target  parameters  estimates.  In  this  work  we  exploit  the  relation  between  the  ambiguity 
function  and  the  Cramer-Rao  lower  bound  (CRLB)  to  calculate  the  bistatic  CRLBs  of  target 
range  and  velocity  and  so  obtaining  a  local  measure  of  the  estimation  accuracy  of  these 
parameters.  We  also  propose  an  algorithm  for  choosing  in  a  multistatic  scenario,  along  the 
trajectory  of  the  tracked  target,  the  pair  transmitter-receiver  with  the  best  asymptotic 
performance  calculated  in  terms  of  CRLB  on  estimation  accuracy. 


Keywords:  CramerRao  lower  bound,  ambiguity  function,  bistatic  radar,  multistatic  radar 


1  INTRODUCTION 


In  active  radar  a  known  waveform  is  transmitted,  and  the  signal  reflected  from  the  target 
of  interest  is  used  to  estimate  its  parameters.  Typically  the  received  signal  is  a  scaled,  delayed 
and  Doppler-shifted  version  of  the  transmitted  signal.  In  monostatic  configuration,  estimation 
of  the  time  delay  and  Doppler  shift  directly  provides  information  on  range  and  velocity  of  the 
target.  This  is  possible  also  in  the  case  of  bistatic  radar  configuration,  even  if  the  relation 
between  measured  or  estimated  time  delay  and  Doppler  frequency,  and  target  distance  and 
velocity  is  not  linear  [Tsa97].  To  measure  the  possible  global  resolution  and  large  error 
properties  of  the  target  parameters  estimates,  the  ambiguity  function  (AF)  is  often  used  both 
in  mono  and  in  multistatic  scenarios. 

In  this  technical  report,  after  calculating  the  bistatic  ambiguity  function  (BAF)  for  a 
burst  of  linear  frequency  modulated  (LFM)  pulses,  we  exploit  the  relation  between  the  AF 
and  the  Cramer-Rao  lower  bound  (CRLB)  to  calculate  the  bistatic  CRLBs  of  target  range  and 
velocity.  The  bistatic  CRLBs  provide  a  local  measure  of  the  estimation  accuracy  of  these 
parameters.  Moreover,  we  compare  monostatic  and  bistatic  CRLBs  as  a  function  of  the 
number  of  integrated  pulses,  target  direction  of  arrival  (DOA),  and  bistatic  baseline  length. 

The  information  gained  through  the  calculation  of  the  bistatic  CRLBs  can  be  used  in  a 
multistatic  radar  system  for  the  choice  of  the  optimum  transmit-receive  pair.  As  known,  a 
multistatic  radar  system  utilizes  multiple  transmitter  and  receiver  sites  to  provide  several 
different  monostatic  and  bistatic  channels  of  observation,  leading  to  an  increase  in  the 
information  on  a  particular  area  of  surveillance.  The  performance  of  each  bistatic  channel 
heavily  depends  upon  the  geometry  of  the  scenario  and  the  position  of  the  target  with  respect 
to  each  receiver  and  transmitter.  In  this  work  we  approach  the  problem  of  optimally  selecting 
the  transmitter-receiver  (TX-RX)  pair  based  upon  the  CRLB  for  the  bistatic  geometry  of  each 
TX-RX  pair.  The  optimal  pair  is  defined  as  that  exhibiting  the  lowest  bistatic  CRLB  for  the 
target  velocity  or  range.  These  results  can  be  used  for  the  dynamical  selection  of  the  TX-RX 
signals  for  the  tracking  of  a  radar  target  moving  along  a  trajectory  in  a  multistatic  scenario. 


2  THE  AMBIGUITY  FUNCTION 

The  bistatic  geometry  is  pictorially  drawn  in  Figure  1.  The  positions  of  the  transmitters 
receivers,  and  target  are  generic.  Considering  an  ordinary  Cartesian  grid,  the  TX  is  located  at 
point  T,  whose  coordinates  are  (xr,yr),  the  RX  is  located  at  point  R  in  (xR9yR)  and  the 
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target  is  located  at  point  B,  whose  coordinates  are  ( x,y ).  The  triangle  formed  by  the 

transmitter,  the  receiver,  and  the  target  is  called  the  bistatic  triangle. 

As  shown  in  Figure  1,  the  sides  of  the  bistatic  triangle  are  RT  ,  Rr  ,  and  L,  where  RT  is 

the  range  from  the  transmitter  to  the  target,  Rr  is  the  range  from  the  receiver  to  the  target  and 

L  is  the  baseline  between  the  transmitter  and  the  receiver.  The  internal  angles  of  the  bistatic 
triangle,  that  without  lack  of  generality  are  assumed  to  be  positive,  are  a ,  fi ,  and  / .  In 
particular,  the  bistatic  angle  ft  is  the  angle  at  the  apex  of  the  bistatic  triangle,  at  the  vertex 
which  represents  the  target.  Assuming  that  the  coordinates  of  the  transmitter,  receiver,  and 
target  are  known,  it  is  possible  to  calculate  all  the  parameters  of  the  bistatic  triangle.  0T  and 

0R  are  the  look  angle  of  the  transmitter  and  the  look  angle  of  the  receiver,  respectively,  they 

are  measured  positive  clockwise  from  the  vector  normal  to  the  baseline  pointing  forward  the 
target. 


Fig.  1  -  Bistatic  geometry. 
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From  Fig.  1,  we  have  that  0T  =  90°  -  a  ,  0R-y- 90°,  /i  =  1 80°  -  a  -  /3  =  0T  -  0R,  and 
Rt 2  =  7^ 2  +  L2  -  27^T  sin  ^ ,  which  gives  the  range  from  transmitter  to  target  RT ,  as  a 
function  of  the  range  from  receiver  to  target  Rr  and  the  look  angle  of  the  receiver  0R . 

In  Figure  1  we  plot  also  the  target  velocity  vector  V  ;  $  is  the  angle  between  the  target 
velocity  vector  and  the  bistatic  bisector,  which  is  measured  in  a  positive  clockwise  direction 
from  the  bisector.  In  particular,  the  bistatic  bisector  is  represented  by  the  vector  BI ,  where  / 
is  the  in-center  of  the  bistatic  triangle,  whose  coordinates  are  (x7,y7 ) . 

The  coordinates  of  the  incenter  can  be  easily  obtained  as 


(•W/) 


L  +  Rr  +  Rr 


(■ x,y)+ 


R, 


L  +  Rr  +  Rj 


(xT ,  y  t')  + 


R, 


L  +  Rr  +  Rj 


(xr^r) 


In  the  bistatic  geometry,  an  important  parameter  is  the  radial  velocity  Va ,  which  is  the 
target  velocity  component  along  the  bistatic  bisector.  From  the  observation  of  Figure  1,  we 


obtain  V=V- Bh 


BI 


\V  cos^ .  The  bistatic  radar  geometry  can  be  completely  specified  in 


terms  of  any  three  of  the  five  parameters,  0T ,  0R ,  Z,  Rr ,  and  Rt.  In  this  technical  report  we 
will  use  0r,L,  and  Rr. 

We  suppose  that  the  radar  transmits  a  sequence  of  linear  frequency  modulated  (LFM) 
pulses  or  chirps.  The  complex  envelope  of  the  transmitted  unitary  power  signal  is  then  given 
by: 


1  N- 1 

u{t)  —  r—  ^  ux  —  nTR ) , 

n=0 

where 


ux{t) 


-±=exp(jxkt2) 


0 


0  <t<T 
elsewhere 


(1) 


(2) 
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N  is  the  number  of  sub-pulses  for  each  transmitted  burst,  Tr  is  the  burst  repetition  time  and  T 
is  the  duration  of  each  pulse,  with  T<Tr/ 2.  Moreover,  kT 2  =BT  is  the  effective  time- 
bandwidth  product  of  the  signal  and  B  is  the  total  frequency  deviation. 

As  known,  the  ambiguity  function  (AF)  represents  the  time  response  of  a  filter  matched 
to  a  nominal  signal,  when  the  signal  is  received  with  a  delay  za  and  a  Doppler  shift  va 

relative  to  the  nominal  values  zH  and  vH  expected  by  the  filter.  Then,  the  AF  definition  is 
[Tsa97] : 


| x{TH’Ta’vH’ v. )|  =  f  «(* - r>* (t-rH) exp (~j2>r(yH  -v„)t)dt 


(3) 


where  u(t)  is  the  complex  envelope  of  the  signal,  za  and  va  are  the  actual  delay  and  Doppler 
frequency  of  the  radar  target  respectively  and  th  and  vH  are  the  hypothesized  delay  and 
frequency.  In  the  monostatic  case  there  is  a  linear  relation  between  za  and  va  and  the  range 
position  Ra  and  radial  velocity  Va  of  the  target,  za=2Ra!c  and  va  = -2Fa/c/c .  Similar 
relations  hold  for  zH  and  vH.  Clearly,  \x{zH,Ta,vH,va)\  is  maximum  for  za-zH  and 

Va=VH- 

Based  upon  the  definition  (3)  we  can  calculate  the  monostatic  AF  for  the  signal  u(t )  in 
(1)  and  (2)  as  [Lev04]: 


N-l)  (  \r-nT 

Z  exP  [JxHN  - 1 + pWr  ]  i  -  - — =-J- 


sin 

[7ff{v-k{z-pTR)){\- 

sin 

\jiv(^N -\ 

\p\)T*] 

7iT{v-k(z-pTR)){\-\i 

r~pTR\/T) 

sin 

for  |r  -  pTR  |  <  T 


0 


elsewhere 


(4) 


where  we  set  z  -zH-za  and  v  =  vH  -ya . 

If  we  limit  the  delay  to  the  mainlobe  area,  namely  to  |r|  <  T  (p=0),  eqs.  (3)  and  (4) 
reduce  to: 
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\x(t,v)\  = 


'  jzf 

T 


sin 

~nT  1- 

H/r)] 

sinf/rWNfT^] 

I 

l 

\t[ 

It) 

AsinfVvT^) 

for  r  <  T 


(5) 


The  AF  exhibits  its  maximum  in  r  =  0  and  v  =  0 .  The  AF  of  eq.  (4)  is  plotted  in  Fig.  2  for 
BT=  20,  Tr=  1  s,  7M).l  s  and  A=8. 


To  refer  eqs.  (4)  and  (5)  to  the  bistatic  geometry  of  Fig.  1  and  to  obtain  the  expression  of 
the  bistatic  AF,  we  must  replace  in  (4)-(5)  the  relations  [Tsa97]: 

/D  n  t^_Rr+\IRr+l2  +  2RrL sin 0R  ^ 

^ H  (Rr  ’  ^R  •>  (^) 

c 

and 
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(7) 


fa 


1 


(Rr,Vcos</>,0r,L)  =  2^Vcos</>  —  + 


Rr  +  L  sin  ^ 


2^Rl+L2+2RRLsm0R 


The  plots  of  eq.  (4)  are  reported  in  Figures  3-7  for  different  values  of  BT  and  TR  in  the  bistatic 
plane  Rr-Vb  with  VB  =  V  cos  $  with  Va=600  m/s,  Ra= 60  Km.  For  low  values  of  BT  the  shape 

of  the  bistatic  ambiguity  function  heavily  depends  on  the  target  angle  0R ,  as  evident 
comparing  Figs. 3-4  to  Fig. 5.  For  some  particular  values  of  0R  as  0R  -  /r/2 ,  the  range 
resolution  is  very  poor  (see  Fig.  5c).  For  high  values  of  BT  the  ambiguity  function  seems  to  be 
less  dependent  on  0R ,  as  apparent  in  Figs.  6  and  7. 

For  completing  the  analysis  in  Figs.  8-10  the  slices  of  AF  are  plotted  for  N=  1,  2  and  32 
showing  that  for  increasing  N  the  range  resolution  improves  but  many  peaks  appear  in  the 
bistatic  AF  shape. 
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Fig.  5a  -  Logarithm  of  the  bistatic  AF,  BT=  20,  7^=  Is,  7M).ls,  N=  8,  0R  -  -0.499^ , 

Fa=600m/s,  i?a=60Km 
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Fig.  6a  -  Log  of  the  bistatic  AF,  BT= 2500,  Jfl=lms,  r=250|us,  N=  8,  0R-n ,  Ffl=600m/s, 

^a=60Km 
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Fig.  6b  -  Log  of  the  bistatic  AF,  BT= 2500,  TR=  lms,  T=250|us,  N=  8,  0R-n ,  Fa=600m/s, 

^a=60Km 
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Fig.  6c  -  Log  of  the  bistatic  AF,  BT= 2500,  Jfl=lms,  r=250|us,  N=  8,  0R-n ,  Fa=600m/s, 

^a=60Km 


Fig.  7a  -  Log  of  the  bistatic  AF,  BT= 2500,  T^lms,  r=250|us,  A=8,  0R  -  -0.499;r , 

Fa=600m/s,  ^a=60Km 
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Fig.  7b  -  Log  of  the  bistatic  AF,  BT= 2500,  Ifl=lms,  r=250|us,  7V=8,  0R  -  -0.499^ , 

L^bOOm/s,  i?a=bOKm 
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Fig.  7c  -  Log  of  the  bistatic  AF,  BT= 2500,  T^lms,  r=250|us,  N=  8,  0R  -  -0.499;r , 

Fa=600m/s,  ^a=60Km 
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Fig.  8  -  Log  of  the  bistatic  AF,  BT= 2500,  7^=1  ms,  r=250jus,  N=  1,  0R  « -;r/2 ,  Ffl=600m/s, 

^a=60Km 
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Fig.  9  -  Log  of  the  bistatic  AF,  BT= 2500,  7)?=  lms,  r=250jus,  77=2,  0R  ^  -/r/2 ,  Fa=600m/s, 

7?a=60Km 
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Fig.10  -  Log  of  the  bistatic  AF,  BT= 2500,  7^=1  ms,  r=250|us,  N=  32,  0R  «  -;r/2 ,  Fa=600m/s, 

i?a=60Km 
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3  BISTATIC  CRAMER-RAO  LOWER  BOUNDS 

Unlike  the  AF  which  provides  information  on  the  global  resolution,  the  CRLB  provides 
a  local  measure  of  estimation  accuracy.  Both  can  be  used  to  asses  the  error  properties  of  the 
estimates  of  the  signal  parameters.  In  [Van71]  the  author  derived  a  relationship  between 
CRLB  and  AF,  which  has  been  successfully  used  in  the  analysis  of  passive  and  active  arrays 
[DogOl].  In  the  monostatic  configuration,  [Van71]  claims  that  for  the  Fisher  Information 
Matrix  (FIM)  the  following  relationship  holds: 


3M(ra,va)  =  -2SNR 


52h(r,^)| 

dr2 

92h(r,vQ| 

dvdr 


d2  |x(r,i/)| 
drdv 

d2  \x (r,v)| 


(8) 


where  SNR  is  the  signal-to-noise  power  ratio  at  the  receiver.  In  Appendix  A  we  report  the 
proof  of  relation  (8).  From  eq.  (8)  the  CRLBs  follow:  CRLB{ra)  =  \JM(Ta’ya)\  I  an^ 

CRLB{va)  =  [3M{Ta,va)\l2. 

In  the  bistatic  configuration  we  should  write  the  AF  in  terms  of  the  bistatic  t(Rr,Or,L ) 
and  v(Rr,Vb,0r,L)  and  derive  it  with  respect  to  the  useful  parameters  Rr  and  Vr.  Then 


J  b(Rr,Vb)  =  -2SNR 


d2\x(RR,VB)\  82\x(Rr,Vb)\ 


8R, 


8Rr8Vb 


82\X(RR,VB)\  82\X(RR,VB)\ 


8Vb8Rr 


8Va 


(9) 


Rn=R„,V„=V„ 


After  some  algebra,  it  is  possible  to  verify  that,  in  the  monostatic  configuration  we  have: 


K(ra,va)  =  ~2SNR 
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then 


CRLB(za)  = 


3 

2ji1T1k1SNR 


r  T\2 


T 

'  i<  J 


1 

TV2  -1 


and 


CRLB(va)  = 


_ 3 _ 

2tt2TrSNR[n2  -l) 


(11) 


(12) 


These  results  are  in  agreement  with  those  obtained  in  [Dog07]. 

For  the  calculation  of  the  CRLBs  in  the  bistatic  domain  we  can  partially  use  the  results 
for  the  monostatic  domain.  Following  the  “chain  rule”  of  the  derivative  (see  Appendix  B  for 
details)  we  can  prove  that 
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From  eq.  (6)-(7)  we  get: 
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dv 

dRR 


f 

—  Vn 
2c 


L2  cos2  0l 


{^Rr  + 1}  +  2RrL  sin  0R  ^ 


3/2  1 


Rr+L  sin  0R 


2  2^7? 2  +  Is  +  2RrL  sin  dh 


(16) 


tv  =  2/’.  jl  |  7?^  +  Z sin dR 

tVn  c  ]j  2  2^R]  +  [}  +2RRLsm0R 


(17) 


07 

5^ 


c 

V 


\ 

Rr+L  sin  0R 
^ Rr  +  Z/2  +  2RrL  sin  0R  ^ 
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(20) 


- 

6^Rr  +  Z/2  +  2  RrL  sin  0R  j  ^ Rr  +  L  sin  0R  ^  ^ Rr  + 12  +  2RrL  sin  ^  ^  +  5  +  L  sin  0R  ^ 

4c  | 

( V2  fcVBL 2  cos2  0R ( R2r  +  L2+  2RrL  sin  0R 

1/2  ~ 
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■3/2 
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8t  _  32t  _  82z  _  82v 

~8Vb~~8rJwb  ~8Vb2~8Vb2 


(22) 


For  L=0,  bistatic  CRLBs  coincide  with  monostatic  CRLBs. 


4  NUMERICAL  ANALYSIS 


The  CRLBs  are  plotted  in  Figs.  11-14  as  a  function  of  the  baseline  length  L,  the  target 
angle  0R ,  and  the  number  of  subpulses  N  in  each  transmitted  burst.  It  is  evident  that,  for  all 
the  parameter  values  we  tested,  the  bistatic  CRLBs  (CRLBb)  are  higher  than  the  monostatic 
CRLBs  (CRLBm).  The  bistatic  CRLBs  get  even  worse  for  target  angles  0R  close  to  -nil, 

where  they  tend  to  infinity.  This  behavior  is  in  agreement  with  the  shape  of  the  bistatic  AF 
plotted  in  Fig.5c,  where  the  AF  exhibit  and  almost  step-wise  shape. 

To  highlight  the  differences  between  monostatic  and  bistatic  domain  we  plot  in  Figs. 
15-22  the  log  of  the  ratio  between  CRLBb  and  CRLBM  as  a  function  of  the  baseline  length  L 
and  the  true  position  of  the  target  Ra.  For  high  values  of  BT  we  observe  that  generally  the 
CRLBb  are  higher  than  the  CRLBm,  particularly  for  the  range  and  when  the  distance  between 
the  receiver  and  the  target  is  orthogonal  to  the  baseline.  When  BT  is  low,  the  CRLBb  show 
some  gain  with  respect  to  the  monostatic  scenario,  even  for  0R=-n/2  as  shown  in  Figs.  19 
and  21. 
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Fig.  11  -  RCRLB  as  a  function  of  bistatic  baseline  length  L,  BT= 2500,  r^=lms,  T=250ps, 
N=  32,  0R  =  n/2 ,  Va=600  m/s,  £a=60Km,  SNR=0dB 
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Fig.  12  -  RCRLB  as  a  function  of  bistatic  baseline  length  L,  BT= 2500,  r^=lms,  r=250|us, 
N=32,  6r=tt  ,  Va=600  m/s,  Ra=60Km,  SA/R=0dB. 
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Fig.  13  -  RCRLB  as  a  function  of  target  angle  0R ,  Z=10Km,  BT= 2500,  7>=lnis,  r=250|us, 

N=32,  Va=600  m/s,  Ra=60Km,  $VR=0dB 
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Fig.  14  -  RCRLB  as  a  function  of  the  number  of  pulses  N,  Z=10Km,  BT=2500 ,  7i=lms, 
T=250ps,  Va=600  m/s,  Ra=60Km,  0R  =  tt/6  ,  SVR=0dB. 
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Fig.  15  -  Log-ratio  between  CRLBb  and  CRLBM  for  the  range,  BT=2500 ,  TR  =lms,  T=250ps, 

7V=32,  eR  =  -71  / 2 ,  Va  =600m/s. 
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Fig.  16  -  Log-ratio  between  CRLBb  and  CRLBM  for  the  velocity,  BT= 2500,  TR  =lms, 
r=250jas,  7V=32,  0R  =  -n) 2 ,  =600m/s. 
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Fig.  18  -  Log-ratio  between  CRLBb  and  CRLBM  for  the  velocity,  BT= 2500,  TR  =lms, 

r=250|us,  N=32,  0R=7r  ,Va  =600m/s 
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Fig.  19  -  Log-ratio  between  CRLBb  and  CRLBM  for  the  range,  BT=  40,  TR  =ls,  7=0. Is,  7V=8, 

dR=  n  ,Va  =600m/s 
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Fig.  20  -  Log-ratio  between  CRLBb  and  CRLBM  for  the  velocity,  BT=  40,  TR  =ls,  T=0.1s, 

7V=8,  0R-n  ,Va  =600m/s 
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Fig.  21  -  Log-ratio  between  CRLBb  and  CRLBM  for  the  range,  BT=  40,  TR  =ls,  7=0. Is,  iV=8, 


0R  =  -nj 2 ,  Va  =600m/s 


-25- 


L  (m) 


x  104 


Ra(m) 


1.8 

1.6 

1.4 

1.2 

1 

0.8 

0.6 

0.4 

0.2 

0 

-0.2 


Fig.  22  -  Log-ratio  between  CRLBb  and  CRLBM  for  the  velocity,  BT=  40,  TR  =ls,  TMfls, 

7V=8,  eR  =  -n  12 ,  Va  =600m/s. 
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5  OPTIMAL  SELECTION  OF  THE  TX-RX  PAIR  IN  A  MULTISTATIC  RADAR 

NETWORK 

The  CRLB  study  carried  out  on  the  bistatic  geometry  can  be  applied  for  the  selection  of 
the  TX-RX  pair  in  a  multistatic  radar  network.  Multistatic  radar  utilises  multiple  transmitters 
and  receivers  to  provide  several  different  monostatic  and  bistatic  channels  of  observation, 
leading  to  an  increase  in  the  information  on  a  particular  area  of  surveillance.  The  information 
gain  obtained  through  this  spatial  diversity  can  give  rise  to  a  number  of  advantages  over  both 
the  individual  monostatic  and  bistatic  cases  in  typical  radar  functions,  such  as  detection, 
parameter  estimation,  tracking  and  identification. 

We  verified  in  previous  sessions  that  the  performance  of  each  bistatic  channel  heavily 
depends  upon  the  geometry  of  the  scenario  and  the  position  of  the  target  with  respect  to  each 
receiver  and  transmitter.  In  this  paragraph  we  approach  the  problem  of  optimally  selecting  the 
transmitter-receiver  (TX-RX)  pair,  based  on  the  information  provided  by  the  CRLB  for  the 
bistatic  geometry  of  each  TX-RX  pair.  The  best  pair  is  defined  as  that  exhibiting  the  lowest 
bistatic  CRLB  for  the  target  velocity  or  range.  These  results  can  be  used  for  the  dynamical 
selection  of  the  TX-RX  signals  for  the  tracking  of  a  radar  target  moving  along  a  trajectory  in  a 
multistatic  scenario. 

In  our  scenario  we  considered  a  map  of  dimension  Lx  -  20 km  and  Ly  -  20 km  and  we 

placed  five  transmitters  and  four  receiver  in  this  area.  In  particular,  we  placed  the  transmitters 
at  coordinates: 

T(1)  =  =  (5km,\5km) 

T{2)  =  (x{2),y{2))  =  (\5km,\5km) 

r(3)  =  (j43),j43))  =  (l0km,10km)  (23) 

J(4)  =  (x{4\y{4)^  =  (5  km,  5km) 
f(5)  =  (-45),j45))  =  (I5km,5km) 

and  the  receivers  at  coordinates: 

R(l)  =  ,  j41})  =  (5km,l0km) 

R{2)  =  (x{2\y{2)^)  =  (I0km,l5km)  (24) 

R(3)  =  (x^3),j43))  =  (\0km,5km) 

R{4)  =  (^4),T^4))  =  (l5km,\0km) 


-27- 


Therefore,  there  are  NT  x  NR  =  5  x  4  =  20  pairs  of  transmitter-receiver  that  we  considered  as 

independent  bistatic  systems.  The  20  resulting  pairs  of  Tx-Rx  are  listed  in  Table  1. 

We  assume  that  each  transmitter  sends  a  burst  of  N  =  8  chirp  pulses  with  a  compression 
ratio  BT  -  2500  and  a  pulse  repetition  interval  PRI  =  10-3  sec .  The  carrier  frequency  of  the 
system  is  fc  =  3  •  1  08/2tt  F/z. 

For  each  point  of  the  analyzed  area  and  for  each  of  the  20  bistatic  systems  we  calculated 
the  root  Cramer-Rao  lower  bound  (RCRLB)  of  the  target  range  and  of  the  target  velocity.  In 
particular,  we  assumed  that,  in  each  point  of  the  analyzed  area,  the  target  has  a  velocity  vector 
aligned  to  the  x  axis  and  with  intensity  of  500  ml  sec. 

As  verified,  the  RCRLB  of  the  target  range  and  of  the  target  velocity  are  a  function  of 
the  range  from  receiver  to  target  Rr  ,  the  baseline  Z,  the  look  angle  of  the  receiver  0R ,  the 

radial  velocity  Va  and  the  signal-to-noise  power-ratio  (SNR).  All  these  parameters  depend  on 

the  configuration  of  the  bistatic  triangle,  that  is,  on  the  coordinates  of  the  target,  the 
transmitter  and  the  receiver.  Bistatic  geometry  also  affects  the  received  echo  power,  because 

the  path  loss  factor  in  this  case  is  (RtRr  )  [SkoOl].  In  particular  the  SNR  can  be  written  as: 


SNR  = 


(35) 


where  SNRC  is  a  constant  parameter.  We  assumed  that  SNRC  =10 dB ,  that  is,  we  assumed 


that  if  both  the  transmitter  and  the  receiver  are  located  in  (0,0)  and  the  target  is  located  in 


(Lx,Ly) 1 ,  then  SNR  =  \0dB. 

Figures  23-42  are  colour  coded  maps  representing  the  RCRLBb  of  the  target  range  and 
of  the  target  velocity  in  each  point  of  the  analyzed  area.  In  particular,  Figure  z'-( a)  represents 
the  RCRLBb  of  the  target  range,  measured  in  dB(m) ,  for  the  zth  bistatic  system;  while  Figure 
z'-(b)  represents  the  RCRLBb  of  the  target  velocity,  in  dB(m  /  sec) ,  for  the  zth  bistatic  system. 

As  apparent  from  the  results,  the  RCRLB  of  each  bistatic  system  is  strongly  related  to 
the  bistatic  geometry.  It  is  clear  that  the  effects  of  geometry  factors  are  more  prominent  as  the 
target  approaches  the  baseline,  that  is,  when  Rr<L  and  the  receiver  look  angle  0R 
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approaches  —nj 2 .  Moreover,  when  the  target  is  on  the  baseline,  the  resulting  delay  is  Lie  and 
the  radial  velocity  is  zero,  and  therefore  resolution  is  totally  lost.  However,  it  is  evident  that 
the  effects  of  the  bistatic  geometry  are  less  noticeable  when  the  distance  to  the  target 
increases;  in  this  case  the  bistatic  system  behaves  more  and  more  as  a  monostatic  system 
Therefore,  the  performances  of  each  bistatic  system  are  strongly  related  to  the  configuration 
of  the  bistatic  triangle,  that  is,  to  the  positions  of  the  transmitter,  the  receiver  and  the  target. 

It  is  clear  that  using  different  transmitting  and  receiving  systems,  the  target  can  be  seen 
by  different  bistatic  configuration;  therefore,  knowing  the  coordinates  of  each  transmitter  and 
each  receiver  of  the  whole  system,  it  is  possible  to  calculate,  for  each  point  of  the  analyzed 
area,  which  is  the  transmitter-receiver  pair  having  the  best  performances,  that  is  the  minimum 
CRLBb. 

Figures  43 -(b)  and  44-(b)  show,  in  a  colour  coded  map,  the  transmitter-receiver  pair 
which  has  the  minimum  RCRLBb  for  each  point  of  the  analyzed  area.  The  colormap  of  these 
figures  in  divided  into  20  colours,  each  of  which  is  associated  with  one  of  the  20  bistatic 
systems  listed  in  Table  1. 

Figures  43 -(a)  and  44-(a)  show  the  minimum  RCRLBb  of  the  target  range  and  the  target 
velocity,  respectively,  that  is  the  value  of  the  RCRLB  which  is  provided  by  the  transmitter- 
receiver  pair  which  has  the  minimum  RCRLB. 

Figures  45-66  show  same  results  as  in  Fig.  23-42  but  obtained  by  choosing  different 
transmitter-receiver  configurations.  In  particular,  the  positions  of  the  5  transmitters  and  the  4 
receivers  have  been  randomly  chosen  inside  the  analyzed  area. 


1  Lx  and  Ly  are  the  projections  of  the  baseline  L  on  the  x  and  y  axis  respectively 
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Table  1:  Analyzed  bistatic  systems. 
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Figure  23-(a)  -  Bistatic  RCRLB  of  the  target  range  [dB(m)].  Pair  1. 


Figure  23-(b):  Bistatic  RCRLB  of  the  target  velocity  [dB(m/sec)].  Pair  1. 
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Figure  24-(a):  Bistatic  RCRLB  of  the  target  range  [dB(m)].  Pair  2. 


Figure  24-(b):  Bistatic  RCRLB  of  the  target  velocity  [dB(m/sec)].  Pair  2. 
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Figure  25-(a):  Bistatic  RCRLB  of  the  target  range  [dB(m)].  Pair  3. 


Figure  25-(b):  Bistatic  RCRLB  of  the  target  velocity  [dBm/sec].  Pair  3. 
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Figure  26-(a):  Bistatic  RCRLB  of  the  target  Range  [dBm].  Pair  4. 


Figure  26-(b):  Bistatic  RCRLB  of  the  target  velocity  [dBm/sec].  Pair  4. 
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Figure  27-(a):  Bistatic  RCRLB  of  the  target  Range  [dBm].  Pair  5. 


Figure  27-(b):  Bistatic  RCRLB  of  the  target  velocity  [dBm/sec].  Pair  5. 
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Figure  28-(a):  Bistatic  RCRLB  of  the  target  Range  [dBm].  Pair  6. 


Figure  28-(b):  Bistatic  RCRLB  of  the  target  velocity  [dBm/sec].  Pair  6. 
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Figure  29-(a):  Bistatic  RCRLB  of  the  target  Range  [dBm].  Pair  7. 


Figure  29-(b):  Bistatic  RCRLB  of  the  target  velocity  [dBm/sec].  Pair  7. 
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Figure  30-(a):  Bistatic  RCRLB  of  the  target  Range  [dBm].  Pair  8. 


Figure  30-(b):  Bistatic  RCRLB  of  the  target  velocity  [dBm/sec].  Pair  8. 
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Figure  31-(a):  Bistatic  RCRLB  of  the  target  Range  [dBm].  Pair  9. 


Figure  31-(b):  Bistatic  RCRLB  of  the  target  velocity  [dBm/sec].  Pair  9. 
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Figure  32-(a):  Bistatic  RCRLB  of  the  target  Range  [dBm].  Pair  10. 


Figure  32-(b):  Bistatic  RCRLB  of  the  target  velocity  [dBm/sec].  Pair  10. 
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Figure  33-(a):  Bistatic  RCRLB  of  the  target  Range  [dBm].  Pair  11. 


Figure  33-(b):  Bistatic  RCRLB  of  the  target  velocity  [dBm/sec].  Pair  11. 
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Figure  34-(a):  Bistatic  RCRLB  of  the  target  Range  [dBm].  Pair  12. 


Figure  34-(b):  Bistatic  RCRLB  of  the  target  velocity  [dBm/sec].  Pair  12. 
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Figure  35-(a):  Bistatic  RCRLB  of  the  target  Range  [dBm].  Pair  13. 


Figure  35-(b):  Bistatic  RCRLB  of  the  target  velocity  [dBm/sec].  Pair  13. 
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Figure  36-(a):  Bistatic  RCRLB  of  the  target  Range  [dBm].  Pair  14. 


Figure  36-(b):  Bistatic  RCRLB  of  the  target  velocity  [dBm/sec].  Pair  14. 
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Figure  37-(a):  Bistatic  RCRLB  of  the  target  Range  [dBm].  Pair  15. 


Figure  37-(b):  Bistatic  RCRLB  of  the  target  velocity  [dBm/sec].  Pair  15. 
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Figure  38-(a):  Bistatic  RCRLB  of  the  target  Range  [dBm].  Pair  16. 


Figure  38-(b):  Bistatic  RCRLB  of  the  target  velocity  [dBm/sec].  Pair  16. 
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Figure  39-(a):  Bistatic  RCRLB  of  the  target  Range  [dBm].  Pair  1 7. 


Figure  39-(b):  Bistatic  RCRLB  of  the  target  velocity  [dBm/sec].  Pair  1 7. 
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Figure  40-(a):  Bistatic  RCRLB  of  the  target  Range  [dBm].  Pair  18. 


Figure  40-(b):  Bistatic  RCRLB  of  the  target  velocity  [dBm/sec].  Pair  18. 
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Figure  41-(a):  Bistatic  RCRLB  of  the  target  Range  [dBm].  Pair  19. 


Figure  41-(b):  Bistatic  RCRLB  of  the  target  velocity  [dBm/sec].  Pair  19. 
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Figure  42-(a):  Bistatic  RCRLB  of  the  target  Range  [dBm].  Pair  20. 


Figure  42-(b):  Bistatic  RCRLB  of  the  target  velocity  [dBm/sec].  Pair  20. 
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Figure  43-(a):  Minimum  RCRLB  of  the  target  range  [dBm]. 


Figure  43-(b):  Optimum  pair  map  for  target  range  estimation. 
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Figure  44-(a):  Minimum  RCRLB  of  the  target  velocity  [dBm/sec]. 


Figure  44-(b):  Optimum  pair  map  for  target  velocity  estimation. 
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Figure  45-(a):  Bistatic  RCRLB  of  the  target  Range  [dBm].  Pair  7,  2nd  configuration 


Figure  45-(b):  Bistatic  RCRLB  of  the  target  velocity  [dBm/sec].  Pair  7,  2nd  configuration 
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Figure  46-(a):  Bistatic  RCRLB  of  the  target  Range  [dBm].  Pair  2,  2nd  configuration 


Figure  46-(b):  Bistatic  RCRLB  of  the  target  velocity  [dBm/sec].  Pair  2,  2nd  configuration. 
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Figure  47-(a):  Bistatic  RCRLB  of  the  target  Range  [dBm].  Pair  3,  2nd  configuration 


Figure  47-(b):  Bistatic  RCRLB  of  the  target  velocity  [dBm/sec].  Pair  3,  2nd  configuration 
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Figure  48-(a):  Bistatic  RCRLB  of  the  target  Range  [dBm].  Pair  4 ,  2nd  configuration 


Figure  48-(b):  Bistatic  RCRLB  of  the  target  velocity  [dBm/sec].  Pair ,  2nd  configuration. 
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Figure  49-(a):  Bistatic  RCRLB  of  the  target  Range  [dBm].  Pair  5.  2nd  configuration 


Figure  49-(b):  Bistatic  RCRLB  of  the  target  velocity  [dBm/sec].  Pair  5.  2nd  configuration 
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Figure  50-(a):  Bistatic  RCRLB  of  the  target  Range  [dBm].  Pair  6.  2nd  configuration 


Figure  50-(b):  Bistatic  RCRLB  of  the  target  velocity  [dBm/sec].  Pair  6.  2nd  configuration 
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Figure  51-(a):  Bistatic  RCRLB  of  the  target  Range  [dBm].  Pair  7.  2nd  configuration 


Figure  51-(b):  Bistatic  RCRLB  of  the  target  velocity  [dBm/sec].  Pair  72nd  configuration. 
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Figure  52-(a):  Bistatic  RCRLB  of  the  target  Range  [dBm].  Pair  8.  2nd  configuration 


Figure  52-(b):  Bistatic  RCRLB  of  the  target  velocity  [dBm/sec].  Pair  8.  2nd  configuration 
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Figure  53-(a):  Bistatic  RCRLB  of  the  target  Range  [dBm].  Pair  9.  2nd  configuration 


Figure  53-(b):  Bistatic  RCRLB  of  the  target  velocity  [dBm/sec].  Pair  9.  2nd  configuration 
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Figure  54-(a):  Bistatic  RCRLB  of  the  target  Range  [dBm].  Pair  10.  2nd  configuration 


Figure  54-(b):  Bistatic  RCRLB  of  the  target  velocity  [dBm/sec].  Pair  10.  2nd  configuration 
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Figure  55-(a):  Bistatic  RCRLB  of  the  target  Range  [dBm].  Pair  11.  2nd  configuration 


Figure  55-(b):  Bistatic  RCRLB  of  the  target  velocity  [dBm/sec].  Pair  11.  2nd  configuration 
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Figure  56-(a):  Bistatic  RCRLB  of  the  target  Range  [dBm].  Pair  12.  2nd  configuration 


Figure  56-(b):  Bistatic  RCRLB  of  the  target  velocity  [dBm/sec].  Pair  12.  2nd  configuration 
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Figure  57-(b):  Bistatic  RCRLB  of  the  target  velocity  [dBm/sec].  Pair  13.  2nd  configuration 


Figure  58-(b):  Bistatic  RCRLB  of  the  target  velocity  [dBm/sec].  Pair  14.  2nd  configuration 
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Figure  59-(a):  Bistatic  RCRLB  of  the  target  Range  [dBm].  Pair  15.  2nd  configuration 


Figure  59-(b):  Bistatic  RCRLB  of  the  target  velocity  [dBm/sec].  Pair  15.  2nd  configuration 
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Figure  60-(a):  Bistatic  RCRLB  of  the  target  Range  [dBm].  Pair  16.  2nd  configuration 


Figure  60-(b):  Bistatic  RCRLB  of  the  target  velocity  [dBm/sec].  Pair  16,  2nd  configuration. 
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Figure  61-(a):  Bistatic  RCRLB  of  the  target  Range  [dBm].  Pair  17.  2nd  configuration 


Figure  61-(b):  Bistatic  RCRLB  of  the  target  velocity  [dBm/sec].  Pair  17.  2nd  configuration 
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Figure  62-(a):  Bistatic  RCRLB  of  the  target  Range  [dBm].  Pair  18.  2nd  configuration 


Figure  62-(b):  Bistatic  RCRLB  of  the  target  velocity  [dBm/sec].  Pair  18.  2nd  configuration 
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Figure  63-(a):  Bistatic  RCRLB  of  the  target  Range  [dBm].  Pair  19.  2nd  configuration 


Figure  63-(b):  Bistatic  RCRLB  of  the  target  velocity  [dBm/sec].  Pair  19.  2nd  configuration 
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Figure  64-(a):  Bistatic  RCRLB  of  the  target  Range  [dBm].  Pair  20.  2nd  configuration 


Figure  64-(b):  Bistatic  RCRLB  of  the  target  velocity  [dBm/sec].  Pair  20.  2nd  configuration 
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Figure  65-(a):  Minimum  RCRLB  of  the  target  Range  [dBm],  2nd  configuration 


Figure  65-(b):  Optimum  pair  map  for  target  range  estimation,  2nd  configuration 
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Figure  66-(a):  Minimum  RCRLB  of  the  target  velocity  [dBm/sec],  2nd  configuration 


Figure  66-(b):  Optimum  pair  map  for  target  velocity  estimation,  2nd  configuration 
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6  FINAL  REMARKS 

In  this  technical  report,  we  calculated  the  bistatic  AF  for  a  burst  of  LFM  pulses.  We 
exploited  the  relation  between  the  AF  and  the  CRLB  to  calculate  the  bistatic  CRLBs  of  target 
range  and  velocity.  The  bistatic  CRLBs  provide  a  local  measure  of  the  estimation  accuracy  of 
these  parameters.  Then,  we  compared  monostatic  and  bistatic  CRLBs  as  a  function  of  the 
number  of  integrated  pulses,  target  direction  of  arrival  (DOA),  and  bistatic  baseline  length. 
Subsequently,  we  used  the  information  gained  through  the  calculation  of  the  bistatic  CRLBs 
for  the  choice  of  the  optimum  transmit-receive  pair  in  a  multistatic  radar  system.  The 
performance  of  each  bistatic  channel  heavily  depends  upon  the  geometry  of  the  scenario  and 
the  position  of  the  target  with  respect  to  each  receiver  and  transmitter.  We  approach  the 
problem  of  optimally  selecting  the  TX-RX  pair  based  upon  the  CRLB  for  the  bistatic 
geometry  of  each  pair.  The  optimal  pair  was  defined  as  that  exhibiting  the  lowest  bistatic 
CRLB  for  the  target  velocity  or  range.  These  results  can  be  used  for  the  dynamical  selection 
of  the  TX-RX  signals  for  the  tracking  of  a  radar  target  moving  along  a  trajectory  in  a 
multistatic  scenario. 

In  this  work  we  assumed  that  only  one  couple  of  TX-RX  is  active,  neglecting  the 
interference  among  the  other  transmitters  on  the  selected  one.  This  ideal  situation  can  be 
achieved  selecting  orthogonal  waveforms  for  all  the  transmitters  and  matching  dynamically 
the  selected  receivers  to  the  waveform  of  the  selected  transmitter. 

Ongoing  research  focuses  on  dynamic  reconfiguration  of  selected  TX-RX  as  a  function 
of  RCRLBs  as  well  as  tracking  accuracy.  In  this  phase  we  will  take  also  into  account  the 
contemporaneous  presence  of  all  the  transmitted  signals. 
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APPENDIX  A  -  Relation  between  CRLB  and  AF 

In  this  appendix  we  report  the  proof  of  eq.  (8).  The  complex  signal  received  by  the  radar 
is 


r(t)  =  s(t.  A)  +  w(t)  0  <t<T  (A.l) 

where  w(t)  is  a  zero  mean  complex  Gaussian  process  and  s(/,A)  is  the  received  signal  in 

absence  of  noise  which  depends  on  the  vector  A  =  [Al  A2  ...  AN  f  composed  by  the 

unknown,  nonrandom  parameters  that  we  want  to  estimate.  We  let  the  observation  interval 
[0,T]  be  long  enough  to  completely  contain  the  pulse. 

We  approach  the  problem  by  making  a  ^-coefficient  approximation  of  r{t) .  Assuming  that 


is  a  orthonormal  basis,  the  AT-coefficient  approximation  of  r(t)  is  given  by 


K 

rK :(0  =  Xw(0 

i= 1 

where 

T 

r,  =  \r{t)(pXt)dt 
0 


Substituting  Equation  (A.l)  into  this  expression,  we  obtain 


where 


and 


r,  =5j.(A)  +  wi 


T 

sj(A)=  ^s(t ,  A)<pj(t)dt 

o 


T 


o 


(A.2) 


(A.3) 


(A.4) 


(A.5) 


(A.6) 


(A.7) 
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Therefore,  the  samples  r.  are  independent,  identically  distributed  (IID),  complex  Gaussian 
variables  with  variance  a*  and  mean  st (A) ,  in  short  notation  ~  CAT  A),  .  The 

probability  density  function  of  r.  is  given  by: 


;--^(A)|2 


PrM’  A) 


7T& 


(A.  8) 


In  order  to  find  the  CRLB  of  the  non-random  parameters  An ,  where  n  =  1,2,..., N ,  the  first 

step  is  to  find  the  likelihood  function.  With  obvious  notation,  the  likelihood  function  can  be 
written  as 


Ai  M4a]  =  ^(OA(^(0,a) 


n 

M 


KG,. 


h-^-(A)|2 


(A.9) 


Now,  if  we  let  K  ^  oo,  A1  [rz(/),  A]  is  not  well  defined.  In  [Van71]  it  is  shown  that  we  can 

divide  a  likelihood  function  by  anything  that  does  not  depend  on  A  and  still  have  a  likelihood 
function.  In  order  to  avoid  the  convergence  problem,  we  divide  (A.9)  by 


PrM\H0(rAt)\Ho)  =  Yl—Te 


7IG. 


(A.  10) 


before  letting  K  — »  go  .  Because  this  function  does  not  depends  on  A ,  it  is  legitimate  to  divide 
by  it  here.  Let  define  the  likelihood  function 


A  Mo.  A] 


A)  \rK  (0>  A] 

1^0  ) 


(A.  11) 


Substituting  into  this  expression,  cancelling  common  terms  and  taking  the  logarithm,  we 
obtain 
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(A.  12) 


In  h[rK{t),  A]  =  — 1rZ[K(A)|2  -291e{^,.(A)}]  = 

^  w  i= 1 

=  — tI  K(A)|2  _ ri"si(A) -  W  (A)J 


Letting  K  —>  go  we  have 


i  T  i  T  i  T 

In  A[r(7),  A]  - - -  J|^(L  A)|2  dt  +  —  Jr* (t)s(t,A)dt  +  —  Jr(7)/(7,  A)^ 


w  0 


w  0 


w  0 


From  [Kay93]  the  Fisher  Information  Matrix  (FIM)  is  given  by 


\E\ 

r82  lnA[r(t),A]| 

l 

l  dA>  J 

E\ 

r d2  In  A[r(7),  A]| 

1 

[  J 

E\ 

d2  In  A[r(7),  A]  j 

ri 

[  8An8A  J 

52lnA[r(0,A] 


8AI8AN 


d2  In  A[r(t),  A] 


8A; 


where  E  {  •}  is  the  Expectation  operator. 

By  differentiating  Equation  (A.  13)  with  respect  to  A  we  obtain: 


ainA[r(<),A]___l_j&^A)s.fcA)i#(  A) 


8A 


S  dAn 


w  0 


8s' (t,  A) 

34. 


dt  + 


1  Tr  ,  3s(?,A)  1  rr  3s*(t,A) 

H  o  Gn  w  o 


8A. 


that  can  be  wrote  as 


(A.  13) 


(A.  14) 


(A.  15) 
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<91n  A[r(7),  A]  1 


A  J[KO  -  s(t’ A)]  A)  d’ + \  J[KO  -s(t.  A)] 

°w  0  8An  J 


a* 


8s*  (t.  A) 
8A 


dt  (A.  16) 


Differentiating  again  with  respect  to  Am  we  obtain: 


g2lnA[KO,A]=_l^U_a/(/,A)a,(t,A)+ 1  jgMl, 

8AndAm  a2  J[  8Am  8A„  1  J]  8An8Am  ] 


-2  Ji  ar  5/t  L  J  ^  ^ 


tv  0 


cLld/L 


dt 


(A.  17) 


Taking  the  expectation,  we  obtain: 


<32  In  A[r(7),  A]1  1  Trds*(t,  A)  &(/,A)  ,  1  Tcds(t,  A)  as^/,  A) 


- i  —  n  =  -^_  I  dt-—  I  dt 

8A8A  al  l  8A_  8A  al  „J  8A_  8A„ 


w  0  m 


w  0  m 


(A.  18) 


where  we  observed  that 


E  { r(t)  -  s(t ,  A)}  =  E  |w(/)}  =  0 


(A.  19) 


Assuming  that 

s(t,  A)  =  \[~Eu(t,  A)  (A.20) 


where  E  is  the  energy  of  s(t,A)  and  u(t.  A)  is  a  unitary  energy  signal,  Equation  (A.  18)  can 
be  written  as 


[ d2  In  A[r(t),  A]] 
1  8An8Am  ] 


=  -2  SNR 


a2 

8An8Am 


(A.  18) 


where  SNR  =  E / cr2  is  the  Signal  to  Noise  Ratio. 
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T 

It  is  interesting  to  observe  that  j]w(/5  A)|2  dt  is  the  Ambiguity  Function  of  u(t,  A)  evaluated 

o 

T 

around  its  maximum.  Denoting  A)  =  A)|2  dt ,  it  is  possible  to  write 

o 


crlb,.  -[•>"]„ 


(A.  19) 


where 


J  1 


g2Z(  A) 

dA ,2 


1 

2  SNR 


52^(A) 

dA2dA, 


d2Z(  A) 

8AldA2 

d\z(  A) 

a4 


52Z(A) 

dA,dAN 


52Z(  A)  92Z(A) 


(A.  20) 
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APPENDIX  B  -  CRLB  derivation:  Chain  rule 

Suppose  we  need  to  derive  with  respect  to  x  and  y  the  function  z  =  f(u,v)  where 
u  =  h  (x,  y )  and  v  =  g  (x,  y ) .  Then  for  the  chain  rule: 


dz  _df  du  df  dv  dz  _df  du  df  dv 
dx  du  dx  dv  dx  ’  dy  du  dy  dv  dy 


(B.l) 


For  the  second  derivatives  the  following  relations  hold 


dLf(du^  d^ f  dv  du  d2f 
-  J  +2 — - - +  J 


dx2  du' 


\dxy 


dudv  dx  dx  dv‘ 


dv 

\dx y 


g/g2»  g2vg/ 
du  dx2  dx2  dv 


(B.2) 


P. 

& 


d^z 
T 


32/ 


du 


'  du '' 

d2f  dv  du  d2f 

'dv2 

dudv  dy  dy  dv2 

df  d2u  d2v  df 
du  dy  dy2  dv 


(B.3) 


d2z  d  z 


d  f  du  du  +  d  f  dv  du  ^df  d  u 


dxdy  dydx  du2  dy  dx  dudv  dy  dx  du  dxdy 

d2  f  dv  du  d2  f  dv  dv  df  d2v 
+  — - - +  — - - + 


dudv  dx  dy  dv2  dy  dx  dv  dxdy 


(B.4) 
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